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Abstract. Kinetic lattice-gas models display fragile glass behavior, in spite of their 
trivial Gibbs-Boltzmann measure. This suggests that the nature of glass transition 
might be, at least in some cases, understood in purely kinetic or dynamical terms. 
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1. Introduction 

A long debated problem in glass physics concerns the nature of the dynamical ergodicity 
breaking and its relation with the existence of an underlying equilibrium phase 
transition [|1], [^]. In mode-coupling theory the glass transition appears as a purely 
dynamic effect due to an instability of the equation governing the time correlation of 
density fluctuations [§, ||. In particular, mean-field disordered models of structural 
glasses show that glassy features are associated to a rugged free energy landscape 
and that the origin of the dynamical transition is the existence of a large number 
of metastable states which trap the system for an infinite time [|5| ||. On the other 
hand, the lifetime of metastable states in finite dimensional short range models is 
finite, since it is always possible to nucleate, by a thermally activated process, a 
droplet of the stable phase. Therefore the dynamical transition appear as an artifact 
of the mean-field approximation, and in real glasses this transition would be just a 
finite-time kinetic effect, at least on time scales much smaller than the lifetime of 
metastable states. Recently, the close connection between the non-trivial structure of 
Gibbs equilibrium states and the appearance of a persistent glassy dynamics has been 
established for a certain class of systems ||. However, since the dynamical universality 
classes are smaller than the static ones || [UJ, and since salient features of glasses 



are essentially of dynamical nature [II], it is important to understand to what extent 
glassy behaviour depends on the details of microscopic kinetics. A generic microscopic 
mechanism leading to slow relaxation phenomena was suggested some time ago by 



Fredrickson and Andersen [F2J. It is based on kinetic rules involving only a selection 



of the possible configuration changes compatibly with the detailed balance and the 
Boltzmann distribution. A kinetic rule can be so effective that there is no need to 
introduce an energetic interaction between the particles. Although they are physically 
motivated (e.g. by the cage effect mechanism) these kinetic models are not intended to 
describe the realistic dynamics of glasses, but rather to show that the glass transition 
could be, at least in principle, a purely kinetic or dynamical phenomenon. Taking 
advantage of this idea we have explored the limit case of a three dimensional lattice 
gas model defined only by short-range kinetic constraints and by a trivial equilibrium 
measure ||13|| . Remarkably, this finite dimensional model exhibits a fragile glass behavior 
unrelated to the existence of a thermodynamic phase transition O, 14, 15] (for another 



case and the related experimental situation, see |16|, [17| and jig] ). It provides a simple 
example of how the distinction between the ideal (static or dynamic) and the laboratory 
(i.e. kinetic) glass transition can be very subtle and elusive. In the following we present 
some numerical results showing that this model reproduces qualitatively some aspects of 
the glassy phenomenology, such as history dependence, irreversibility effects, power-law 
approach to the asymptotic state, and simple aging behavior. Some related works on 
constrained lattice-gas models are 0, |(], [21], |22| . 
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2. The model 



Our starting point is a kinetic lattice-gas model introduced by Kob and Andersen [13[ . 
The system consists of iV particles in a cubic lattice of size L 3 , with periodic boundary 
conditions. There can be at most one particle per site. Apart from this hard-core 
constraint there are no other static interactions among the particles. At each time step 
a particle and one of its neighbouring sites are chosen at random. The particle moves if 
the three following conditions are all met: 

(i) the neighbouring site is empty; 

(ii) the particle has less than m nearest neighbours; 

(iii) the particle will have less than m nearest neighbours after it has moved. 

The rule is symmetric in time, detailed balance is satisfied and the allowed configurations 
have the same statistical weight in equilibrium. Significant results are obtained when the 
value of m is set to 4. With this simple definition one can proceed to study the dynamical 
behavior of the model at equilibrium. One observes that the dynamics becomes slower 
and slower as the particle density p increases; in particular, the diffusion coefficient of 
the particles, D, vanishes as the density p approaches the critical value p c ~ 0.88, with 
a power law 

D(p)~(p c -p)t>, (1) 

with an exponent ~ 3.1 [T3|] . Since we are interested in the dynamical approach 
to the putative equilibrium state we allow the system to exchange particles with a 
reservoir characterized by a chemical potential p. Therefore, we alternate the ordinary 
diffusion sweeps with sweeps of creation/destruction of particles on a single layer with 
the following Monte Carlo rule: we randomly choose a site on the layer; if it is empty, 
we add a new particle; otherwise we remove the old particle with probability e _/3A \ The 
number of particles is no longer fixed and the external control parameter is 1/p, which 
plays the role of the temperature. The equilibrium equation of state p = p eq (p) is then 
trivially calculated. There is therefore a critical value p c of p defined by p eq (p c ) = p c 
corresponding to the ideal glass transition of the model. In this way we can prepare 
the system in a non equilibrium state by a process analogous to a quench, which is 
represented by a jump in 1/p from above to below 1/ p c . Or, we can let 1/p decrease 
or increase smoothly like in cooling or heating experiments. The situation becomes 
analogous to the canonic case in which one controls the temperature, and the energy 
endeavors to reach its equilibrium value. 

3. Thermodynamics 

Before to study the non-equilibrium regime let us consider the static properties. The 
point is relevant for the question of whether the possible ideal glass transition is 
purely dynamical or is a consequence of an equilibrium transition of some sort. The 



Kinetic glass transition 



4 



Hamiltonian of the model is 

N 

H=- fJl J2 n ii (2) 

i=l 

where = 0, 1 are occupation site variables and /x is the chemical potential. The 
corresponding partition function, for a system of volume V = L 3 , 

Z=(l+e^) V , (3) 

would describe correctly the thermodynamics of the system provided that the 
measure of configurations made inaccessible by the kinetic constraints vanishes in the 
thermodynamic limit. It is possible to convince oneself that the kinetic rules, which 
satisfy detailed balance, allow an initially empty lattice to be progressively filled in, 
leaving only 0(1/ L) empty sites per unit volume. Indeed, it is always possible to find 
a path connecting almost any two allowed configurations, if necessary by letting the 
particles escape one by one by the way they got in. Therefore the Markov process 
generated by the dynamical evolution rule is irreducible on the full manifold of particles 
configurations and the static properties of the model are described by the (Q). In 
particular the state equation and the entropy are respectively given by: 

p=l/(l + e~^), (4) 
S= -plogp-(l-p)log(l-p) . (5) 

Since the static properties of the system are regular as a function of the density or 
the chemical potential, the possible ideal glass transition should appear as a purely 
dynamical effect. The critical value of /i and S corresponding to the threshold density 
p c , can be estimated from the previous equations and they are given by 

/i c ~ 2.0 , S c ~ 0.36 . (6) 



4. History dependence 

A first insight into the nature of the relaxational processes can be gained by studying 
the behavior of one-time observables (energy, specific volume etc.) in a slow annealing 
procedure. We consider a compression experiment in which the inverse chemical 
potential of the reservoir, is slowly decreased at fixed rate from a value 

corresponding to a low density equilibrium configuration up to zero. The simulation 
results presented in the following refer to a system of size 20 3 . In Figure [I] the numerical 
results of the specific volume v = 1/p vs. are compared, for several annealing 

rates, r, with the equilibrium state equation of the system (the smooth curve). In close 
resemblance with the behavior of real glasses these curves exhibits the characteristic 
annealing dependence of one-time observables: after a certain value of the inverse 
chemical potential, 1/7% (r), is reached the dynamics become so sluggish that the system 
is no longer able to follow the annealing procedure; the faster the compression, the 
sooner the system falls off equilibrium. The limit value of v reaches a plateau that 
depends on the compression rate and never seems to cross the critical value v c = l/p c 
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Figure 1. Compression experiment. The smooth curve is the equilibrium state 
equation v — v(l/fi); while the dashed line is the critical value of specific volume 
v c ~ 1.13. The compression rates are (from top to bottom) r — 1/3 ■ 1CP 4 , 10~ 5 , 1CP 6 
in unit of [/i • MC sweep] ~~ 1 . Inset: the same experiment in a system without kinetic 
constraints; here the rate is r = 10~ 4 . 



(the horizontal dashed line). In the inset of fig. |I| we also show the same plot for a 
compression experiment, but this time by removing the dynamical constraints. We see 
that the ordinary lattice gas has no problem in equilibrating at each value of chemical 
potential /z; therefore our "experimental" setup (the way in which the particles reservoir 
and its connection with the system is realized) provide a suitable representation of the 
equilibrium properties of the model. 

Kauzmann's paradox. Once we obtained the experimental equation of state, we can 
evaluate the entropy variation of the reservoir by numerical integration, which is given 
by: 

Situ) = SQh) - J* ^dfi . (7) 

This "calorimetric entropy" in presence of irreversible effects will be different from the 
thermodynamical entropy S eq . Indeed, in fig. |2] we see that when the relaxation time 
exceeds the inverse of the annealing rate the numerical data remain consistently above 
the equilibrium curve S = S C q(fi). If one were given only the dynamical data of figs. [L] 
and one would feel tempted to extrapolate the equilibrium specific volume and entropy 
to lower values of Then, given that both these quantities are bounded, one could 
conclude that the Kauzmann "temperature" , defined here as 
1 1 
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Figure 2. Calorimetric entropy as obtained by integration of compression experiment 
data. The smooth curve is the equilibrium entropy S = S(l/fi) while the dashed line 
is the critical value of entropy S c ~ 0.365. The compression rates are (from top to 
bottom) r = 1/3 • lO -4 , 1CT 5 , 1CT 6 . 



is different than zero and therefore that there has to be a static transition. This is 
the usual argument, known as Kauzmann's paradox, according to which the glassy 
state is related to the existence of a thermodynamic phase transition. Of course, here 
there is no such static transition: in this simple case we have access to the whole 
equilibrium curves, which are perfectly analytical though they change concavity rather 
sharply. Irreversibility effects are also evident when we let perform a cycle: in this 
case the specific volume appear to follow a hysteresis loop whose area decreases as the 
compression speed decreases (fig. |3[). 

5. Structural relaxation 

We now turn to the behaviour of the system after a sudden quench to a subcritical 
value 1/p < l/p c - In order to allow the system to reach more rapidly the asymptotic 
regime we perform a "gentle" quench i.e. starting from a configuration with density 0.75 
corresponding to a chemical potential closer to p c . Figure f| shows the time relaxation of 
particles density after a subcritical quench at 1/p = 1/2.2. We see that p never exceeds 
the threshold p c , but rather approaches it like a power law in time: 

Pc -p(t)~r z } (9) 

where t is the time elapsed after the quench and where the exponent z ~ 0.3. Therefore 
the diffusion coefficient D of particles after a subcritical quench vanishes as 

D(t) ~ r c , (10) 
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Figure 4. Time relaxation of the excess density after a subcritical quench, Sp = 
(Pc — p)/(Pc — Po)- The straight line is Sp ~ t~ z with z ~ 0.3. 
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with the exponent £ = z<j> quite close to one. This is closely related to the aging behavior 
observed in the two-time mean-squared displacement of particles, B(t,s). Indeed, in a 
simple minded approach such a quantity would be given by 

B(t,s) = J*cLtD{t) , (11) 

from which follows the simple logarithmic aging 

B(t,s)~ log(t/s), (12) 

in good agreement with the numerical results and the analytical solution of the 
associated singular diffusion model [ 15| . Since the size of the system considered 



here is finite, equilibrium will eventually be reached, (since almost any two allowed 
configurations can be connected by a path of allowed moves), but with times which 
grow fast as L — > oo. 



Activated processes. It is interesting to investigate the role of activated hopping 
processes in the low-temperature phase of glassy systems since they are responsible 
of restoring the ergodicity broken at the glass transition, and it is important to know 
the characteristic time scale on which this equilibration process takes place. As pointed 
out in ref. [|C3|, the activation processes can be simply implemented in this model by 
allowing the violation of the kinetic constraint with a given probability p. Figure [| 
shows an example of a v vs. plot in a compression experiment at fixed annealing 
rate for several values of the activation probabilities, p. As expected, in this case the 
system become able to cross the threshold and, after a certain value of p, p*, it follows 
the full equilibrium curve; we can also see that, for p below p*, the dependence of the 
relaxation time from 1/ fi is not affected by p, since the annealing curves depart from the 
equilibrium one approximately at the same point. The relation between the activation 
probability p* and the equilibration time is better characterized by looking at behaviour 
of the density after a sudden quench. In fig. |6| the relaxation curves in presence of 
activated processes are compared with that one obtained previously for p = (fig. f|). 
If we conventionally defined the ergodicity time, r erg (p), as the time at which the curves 
with p y£ depart from that one at p — 0, it appears that this characteristic time follows 
a power-law, T erg (p) ~ p~ a , with an exponent a ~ 1. A similar result was obtained by 
Castellano and Franz (unpublished). This seems to provide a further evidence of the 
existence of a purely dynamical glass transition in this model. 



6. Conclusions 



To summarize we have shown that three dimensional lattice-gas models defined by short 
range kinetic constraint and trivial equilibrium Boltzmann-Gibbs measure display many 
features of the fragile glass behaviour. Glassy phenomena may therefore have a purely 
dynamical or kinetic origin unrelated to an underlying equilibrium phase transition, 
even in finite dimension and in absence of metastable states. If this kinetic model 
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Figure 5. Compression experiment in presence of activated hopping processes for a 
fixed annealing rate, r = 10~ 5 and different activation probabilities p = I0~ k with 
k = 1, 2, 3, 4. 



Q. 0.9 r 

0.88 1 
0.86 

0.84 r 

0.82 1 

0.8 

0.78 h 



0.76 L-^w#M 



1 








Figure 6. Density relaxation in presence of activated hopping processes: p = 10~ k , 
k = 1, 2, 3, 4, 5 (from left to right); the continuous line represents the data with p = 0. 
The subcritical value of the quench, = 1/2.2, corresponds to an equilibrium density 
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exhibits a true dynamical transition it would provide a microscopic realization, in finite 
dimension, of the mechanism invoked by the ideal mode-coupling theory for the glass 
transition. Of course, it is hard to establish from numerical simulations the existence 
of such a dynamical transition. Indeed, a comparison with the backbone percolation 
problem show that the linear size-dependence of the critical threshold cannot be faster 
than 



l-p c (L)~l/log(logL), (13) 

therefore even if lim^oo p c (L) = 1 (i.e. there is no ideal dynamical glass transition), the 
length-scale over which such a value would be observable is not experimentally accessible. 
In this respect, the emergence in purely kinetic models of a well defined macroscopic 
effective temperature [E3J associated with the violation of the fluctuation-dissipation 



theorem appears quite surprising [23J. Indeed, given the non-holonomic nature of 
kinetic constraints and the trivial Hamiltonian of the model, it would be interesting 
to understand whether a statistical mechanics approach based on the calculation of 
some restricted partition function is able to predict the features of the glassy phase 
and in particular the value of the so called fluctuation-dissipation ratio. A first step in 
this direction would consist in to define a kinetic analogue of the metastable states by 
considering, for example, as metastable those system configurations where all particles 
are blocked by the kinetic constraint; and then find a way to count them. 
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